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Abstract. This paper deals with the analytic derivation of a complete equation of 
state (EOS) for solids from shock- wave data in the range of pressures and temperatures, 
attained by detonation of chemical explosives. It is assumed that stresses behind the 
shock wave are isotropic, and the compressed materials behind the shock front are 
in thermodynamic equilibrium. A single Hugoniot curve determined from shock-wave 
experiments does not provide enough thermodynamic information to specify an EOS. 
The assumptions which, along with the shock-wave data, are used in the present work 
to determine a complete EOS are the Griineisen assumption, the volume dependence 
of the Griineisen parameter, and the Debye model for the specific heat. Consequently, 
the thermodynamic functions of the system under consideration are represented as a 
superposition of a cold (0 K) term corresponding to the cold compressed state and a 
second one corresponding to the thermal vibrations of atoms in the crystal lattice about 
their mean positions. A differential equation is derived for the cold isotherm of the 
internal energy, using the Mie-Griineisen model and the shock Hugoniot as a reference 
curve. An analytic solution to the equation is obtained. Based on this solution, the 
complete, the caloric, the thermal and the incomplete EOS are determined. Also, an 
equation is derived for the shock temperature along the Hugoniot curve. 
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Introduction 

The study of shock-wave propagation in solids has added significantly to our 
understanding of physical processes which take place at high pressures, high 
temperatures and very short times. One particular aspect of these studies is to determine 
the equations of state (EOS) for different materials from shock- wave data. Shock- wave 
experiments, aimed at EOS investigations, render it possible to extend the range of 
pressure-volume data beyond the region that can be reached with conventional static 
pressure experiments and to use much bigger specimens of the investigated material. 

Equations of state of matter, obtained both theoretically and experimentally are of 
immense current importance in the basic and the applied sciences. They describe the 
dependence of thermodynamic properties of matter on its microscopic internal structure 
and have a variety of applications. These include studies of the state of matter in 
the Earth's and planetary interiors, shock-induced chemical reactions and polymorphic 
phase transformations, hypervelocity impact, powder compaction or other applications 
as well as their computer simulation based on hydrodynamic models. In principle, all 
the thermodynamic properties of a system can be predicted once precise knowledge of 
the complete EOS becomes known. But a single Hugoniot curve determined from shock- 
wave experiments does not provide enough thermodynamic information to specify an 
EOS. The assumptions which, along with shock- wave data, are used in this work to 
determine a complete EOS are the Griineisen assumption (1), the volume dependence 
of the Griineisen parameter (2), and the Debye model for the specific heat (3). Also, 
it is assumed that stresses behind the shock wave are isotropic, and the compressed 
materials behind the shock front are in thermodynamic equilibrium. 

Therefore the objective of the work presented herein is to obtain a complete EOS 
for solids from shock-wave data using the Debye - Mie-Griineisen model and the shock 
Hugoniot as a reference curve, and to establish the thermodynamic properties of a system 
with such an EOS. Also, an equation for the shock temperature along the Hugoniot curve 
will be derived. 

1. General formulation 

Formally, the equation of state is a functional relationship among the thermodynamic 
variables for a system in equilibrium. One of the most frequently used thermodynamic 
relations is the thermal EOS, or just the EOS. For a simple thermodynamic system it 
has the form 

f(P,V,T) = 0. (1.1) 

It represents a surface in the PVT space. The PV isotherm (T = const.), isentrope (S 
= const.) and the shock Hugoniot (the locus of all possible states that can be reached 
by using a single shock from a given initial state) are particular curves on this surface. 
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To obtain the form of /, it is convenient to derive the Helmholtz free energy of the 
system 

F = E-TS (1.2) 
and compute P as the volume derivative 

p = -(dF/dV) T . (1.3) 

It is clear that F depends upon the microscopic structure of the solid under 
consideration, which would vary as a function of volume and temperature. At different 
temperatures and densities, the corresponding region of matter will be dominated by 
different interactions. 

In view of the above, we may write the free energy as a superposition of terms 
appropriate to various physical phenomena (Griineisen's assumption |jlj|) 

F(V, T) = F C (V) + F mb (V, T) + F e (V, T) (1.4) 

Here F c is the configurational free energy at OK. F v n, corresponds to the contribution of 
the zero and thermal vibrations of the ions of the crystal lattice. F e is the conduction 
electron thermal excitations contribution. 

Expressions for the internal energy, pressure and other thermodynamic functions 
are similar to Eq. (|1.4j) : 

P = P C + P mb + P e , E = E C + E mb + E e . (1.5) 

The contribution of the conduction electrons thermal excitations manifests itself at 
temperatures > 10 4 K and pressures > 10 2 GPa 0. At lower temperatures and pressures 
it may be neglected and Eqs.( |1.4j ) and ( |1.5| ) take the form 

F = F c + F vib , P = P C + P vib , E = E C + E mb . (1.6) 

The vibrational terms are sums of the quantum zero-degree and the thermal vibrations 
of the crystal lattice. Then Eqs.( |1.6|) may be written as 

F{V, T) = F C (V) + F (V) + F lT (V, T), 

E(V,T) =E c (V) + E (V) + E lT (V,T), (1.7) 
P(V, T) = P C {V) + P (V) + Pit{V, T). 

Fq, Eq and Pq are the quantum zero-degree contributions. They are functions of 
volume only, like the configurational terms. Fit, Fit and Pit are the lattice thermal 
contributions. They depend on volume and temperature. 

The pressures and temperatures attained by detonation of chemical explosives in 
contact with the specimen or explosively accelerated flyer plate methods do not exceed 
10 4 K and 10 2 GPa. These techniques are the most accessible to the general scientific 
community and because of that the most frequently used ones. The behaviour of a 
solid within this range of pressures and temperatures is determined by the lattice 
thermal vibrations. That is why it is often called the phonon region. It looks 
narrow on the background of the highest temperatures and pressures, investigated by 
contemporary physics. But in addition to being easily accessible, it covers a wide variety 
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of processes and phenomena of interest to science and technology. Those mentioned in 
the introduction are to name but few. Consequently, we shall aim our considerations at 
the phonon region and seek the explicit form of Eqs.( [L.7|) . 

From Eq.( |1.2| ) and the fundamental thermodynamic identity 

dE = -PdV + TdS 

it follows 

E c = F c and P c = -dFjdV = -dE c /dV. (1.8) 

The specific form of Eqs.( |1.8| ) will be derived from the Mie-Griineisen equation using 
the shock Hugoniot as a reference curve. 

The calculation of Helmholtz free energy components F (V) and Fit (V, T) will be 
performed within the framework of the Debye model for the specific heats. 



2. Calculation of the zero-degree Kelvin isotherm from shock-wave data 

The shock-wave methods for deriving the cold isotherm are based on the measurement of 
the Hugoniot curve. The quantities directly measured are the kinematic parameters of 
the shock wave - the shock front velocity u s and the particle velocity in the compressed 
region u p . The relation u s (u p ) is the shock Hugoniot. Most substances in the absence 
of phase transitions have a linear shock Hugoniot || 

u s = c + sup (2.1) 

in a wide range of particle velocities. The Hugoniot intercept, cq , and the slope, s, are 
determined from the data by the method of least squares. If rigidity effects and possible 
low pressure phase changes are neglected, the intercept should correspond to the velocity 

1 /2 

of an infinitesimal pressure pulse, or the bulk sound speed, cq =[{dP/ dp)s\ (at P — 
0). Since the slope is linearly related to the pressure derivative of the adiabatic bulk 
modulus, (dB/dP)s, a linear u s -Up Hugoniot then reflects a nearly linear dependence of 
Bs on the pressure. 

The transition from kinematic (u s ,u p ) to thermodynamic (P, V,E) variables is done 
using the laws of conservation of mass, momentum and energy across the shock front [fl 

Pou s = p(u s - u p ), (2.2) 
P - P = pou s u p , (2.3) 

.2" 

p u s = P u p . (2.4) 



(E-E )-^ 



Here E, P and p are the specific internal energy, the pressure and the density behind the 
shock front, and Eq, Pq, po are the values of these quantities ahead of the shock front. 
From Eqs. (|2.2|) and (|2.3|) we obtain 

u s = V [(P-P )/(V -V)} 1/2 , (2.5) 
Up=[(P-P )(V -V)] 1/2 . (2.6) 
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Solving Eqs. fl2.5|) and ( j2.6| ) for P gives 

P = Po + ^u s u p , V = V ChJZlL) . (2.7) 

Equations ( |2.7j ) give the relationship between the kinematic variables u s , and the 
thermodynamic variables P and V. 

Substituting with the right-hand sides of Eqs. fl2.5|) -( [2.6D in Eq.( fOl) gives the 
Hugoniot equation of energy 

E - E = (P + P )(Vo - V)/2. (2.8) 

that defines all states on the (E-P-V) surface that can be reached from an initial state 
(Eq, P , Vo) by a single shock. 

If the linear u s -u p relation holds, the Rankin- Hugoniot equations (|2.7Q and (|2.8|) 
can be used to express pressure and energy as functions of volume along the Hugoniot 
by the following convenient analytic expressions (with Pq and Eq taken to be zero at 
ambient conditions) 

Ph = -f^, (2.9) 
(1 - es) 

Eh = ^ 2 , (2.10) 
2 (1 — es) 

s = 1-V/Vq = 1- P o/p. (2.11) 

Equations ( |2.9D -( PTTD ) and the values of Cq and s summarize all the experimental 
thermodynamic information which is available from shock-wave measurements. 

To compute the cold compression curve we shall employ the Mie-Griineisen equation 
in the form 

P(V, T) - P(V, 0) = ^p- [E(V, T) - E{V, 0)] , (2.12) 

where j(V) is the Griineisen parameter. Equations ( |2.9j , j2.10p give the pressure Ph and 
the specific internal energy Eg on the shock Hugoniot. Since it connects equilibrium 
thermodynamic states we can write 

P H (V) - P(V, 0) = ?Qp- [E H (V) - E(V } 0)] . (2.13) 

We shall assume that 7 is a function only of volume and, specifically, that the 
product p7 is constant. Experimental work on a number of materials 0, as well as 
theoretical considerations || , indicate this to be an adequate approximation. With this 
assumption, "f(V) is given by 

l{V)/V = 70/V0, (2.14) 

where 70 is the thermodynamic value at standard conditions given by 70 = 3acl/C p , 
a is the thermal expansion coefficient, cq the sound speed, and C p the specific heat at 
constant pressure. 
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At T = OK F = E and from Eqs.(|L8D it follows 

P C {V) = P(V, 0) = -dE c (V)/dV. (2.15) 

If we substitute with Eqs.(|2~9|),(pJ0|),(pJ5D in the left-hand side of Eg. (gig) , we 
obtain a differential equation for the zero Kelvin isotherm of the specific internal energy 
E c (e)=E(e,0): 



dE c (e) 
de 



- loEc(e) 



cge(l-7o/2) 
(1 - se) 2 



(2.16) 



E c {e = 0) = E { 



oo- 



This is a first order linear inhomogeneous ordinary differential equation. The functional 
form of E c (e) may be obtained by integrating Eq. (|2.16 ). Unfortunately, the solution 
cannot be given by combinations of elementary functions. It is possible, however, that 
it can be approximated by a power series in e within high precision. In this case the 
cold compression curve E c (e) can be written as the following power series in e 

E c {e)=Y J E, l e\ (2.17) 

i 

where e is the dimensionless volume given by Eq. Q2.11| ). Eqo is obtained from the 
assumed reference state and the remaining coefficients to the fifth order in e are given 
by the expressions: 



E m = 70^00, E 02 = - (eg + 7o^oo) , E 03 = - (4sq 



7o-^oo I ■ 



Eqa = 

24 
1 



2( 7o -9s)^s-7^ 



00 



(2.18) 



E 



05 



120 



2 (7 2 + 97oS - 48s 2 ) c 2 s - 7o 5 ^oo 



From Eqs.( ^.l5| ) and (2.17) we can derive an expression for the zero Kelvin isotherm of 
the pressure: 



(2.19) 



The coefficients Poi and Eoi are not independent. A relation exists between them which 
is given by 

P Qi = (i + l)E 0l+l /V . (2.20) 



3. Contribution of the quantum zero Kelvin and lattice thermal vibrations 
to the EOS 

In accordance with the quasiharmonic approximation the vibrational energy levels of 
a crystal lattice with N atoms may be treated as the energy levels of a system of 3./V 
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independent linear harmonic oscillators. For the logarithm of the partition function of 
such a system we can write 

\nZ = 2Jn2i = 2^ In- . . 3.1 

~t ^i 1 _ ex P {-hvi/kT) 

where ^ is the partition function of the i-th oscillator. 

In the Debye model J7] the crystal lattice is replaced by an isotropic elastic medium 
with a continuous dispersion law function. Further it is assumed that the magnitude 
of frequency in such a system does not exceed a certain boundary value, the Debye 
frequency vd , chosen so that the number of the independent lattice oscillations is 
equal to the total number 3iV of the lattice degrees of freedom. It follows from this 
approximation that we can replace the summation in Eq. (|3.1| ) by integration. Let us 
omit the index % of the oscillations. Then 
j exp {—hv/2kT) 
11 1 - exp (-hv/kT) 

is Inz for all oscillators in the frequency range between v and v+dv. The number of 
these oscillators is equal to q(v)dv, where 

g{y) = 9Nu 2 /u 3 D (3.2) 

is the density of distribution of oscillators among frequencies |7|]. We introduce the 
Debye temperature Op = hvo/k, which depends on volume and is specific for each 
substance, and the new integration variable x = hvjkT. Then Eq.( |3.lD takes the form 

9.w„ l r 9 ° /T 

(Od/TY 

Integration by parts of Eq. ( |3.3| ) gives 

InZ = - 9 -N (6 D /T) - 3N\n [1 - exp (-6 D /T)} + /T J-^. 

The last equation can be written in a more compact form if we introduce the function 



9 1 r o /t 

lnZ = - N (e D /T) - 9N— — 3 ^ x 2 In [1 - exp (-*)] dx. (3.3) 



. 3 r z x 3 dx 

D y z > = i / ~x — 1~' 

z 6 Jo e x — 1 

which is known as the Debye function and z = Op/T. Then the logarithm of the 
partition function takes the form 

In Z = --JV (6 D /T) - 3Nln [1 - exp {—9 D /T)\ + ND (6 D /T) . (3.4) 



From Eq.( |3.4j) and the thermodynamic identities 

F = -kT In Z, E = F~ T(dF/dT) Vt 
p = -(dF/dV) T , C v = (dE/dT) v 

we can derive expressions for the specific free and internal energy, the pressure and the 
specific heat: 

F mb (T,V) = lNk6 D + NkT{3\n[l-exp(-6 D /T)]-D(6 D /T)}, 
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E mb (T,V) = !Nk9 D + 3NkTD(9 D /T), 

C v = 3Nk[D(6 D /T)-(e D /T)D'(e D /T)}, 

where (9/8)M$d is the contribution of the quantum zero oscillations in terms of the 
Debye model, and D' is the derivative of D with respect to 0e>/ T. 

In order to derive P V i b it is necessary to compute the volume derivative of In Z. The 
partition function does not depend explicitly on volume but through the dependence of 
the energy levels on it. It is not possible to obtain this relation from first principles. This 
imposes the use of various approximations. In the case of a solid the volume dependence 
of the energy levels is reduced to volume dependence of the lattice frequencies. This 
relation cannot be obtained from first principles as well. According to Griineisen 

dlnvj/dlnV = -7 (j = 1, 2, 3N), (3.5) 

where Vj(V) are lattice frequencies, and 7 is a positive quantity, one and the same for 
all the 3N normal modes of the crystal lattice. It is assumed that these frequencies, 
and hence, 7, do not depend on temperature, but on volume only. This is sometimes 
referred to as the quasiharmonic approximation. 

Equation ( |3.5| ) is the statistical definition of the Griineisen parameter. It holds for 
any frequency, hence it may be written for the Debye frequency i/p 

d\n.V£)/ d\nV = —7 
and since vd and 9d are proportional, it follows 

d\n6 D /d\nV = -7 or d\n6 D /dV = -7/V. 
Accordingly, we obtain for P vib 

Pvib = lh (V) /V] Nk6 D + 3 [7 (V) /V] NkTD(6 D /T) 

o 

or 

Pvib=b(V)/V]E vib , 

which in terms of the assumed approximation 7/V = 70/V0 = const, takes the form 
Pvib = (lo/V )E vib . 
These results render it possible to write Eqs. (|1.7| ) in an explicit form: 

F(T, e) = J2 EoiS* + lNk6 D + NkT {3 In [1 - exp {-9 D /T)\ - D(6 D /T)} , 

i 8 

E(T, e) = V E 0i e l + ^NkOn + 3NkTD(6 D /T), 
i 8 

P(T,e) = Y.P^ 1 + %o/V )Nke D + 3( l0 /V )NkTD(9 D /T), 
where the coefficients Poi are given by Eq. ( p.20|) . 
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These are the complete, the caloric and the thermal EOS for solids. They are 
defined over domains of the (P,e) and (E,e) planes formed by the set of experimental 
values of e, P and E. 

The integration constant E 00 can be easily obtained from the caloric EOS. Let us 
write it in explicit form: 

E(T, e) = E 00 + E 01 e + E m e 2 + ... + %Nk6 D + 3NkTD(9 D /T). 

o 

Since ambient conditions are assumed as a reference state, e.g. E( T=300K,e=0) = 0, 
we have 

E 00 = ~-Nkd D - 900NkD(e D /300). 
8 

Knowing the complete EOS, we are able in principle to derive all the thermodynamic 
properties of the system under consideration by simple differentiation. 

To calculate the temperature on the shock Hugoniot we may use either the caloric 
or the thermal EOS. Let us write the thermal EOS in the form T = f(T) and replace 
P{T,e) with Pff(e). We obtain 



T H {e) = (Vo/3 7 o) 



P H (e) - ]TP ^ ~ %o/V )Nk9 D 



/NkD{9 D /T H ). 



This equation may be solved either graphically or by the method of successive 
approximations. 

The simultaneous solution of the caloric and the thermal EOS in order to eliminate 
T defines the incomplete EOS 

P(E, e) = (jo/V )E + + / V °\ E oi+i ~ (lo/V )E 0i } e\ 

i 

This is the constitutive relation necessary to supplement the system of equations 
of fluid dynamics describing adiabatic compressible inviscid flow. 



4. Conclusions 



The assumptions that the free energy of a solid may be represented as a superposition 
of terms, corresponding to the cold compressed state, the quantum zero vibrations and 
the thermal vibrations of the crystal lattice, that the ratio / ~f(V)/V remains constant, 
and the Debye model for the specific heat, along with the experimental shock Hugoniot 
allow the derivation of its thermodynamic properties in explicit form. In the beginning 
a differential equation for the zero Kelvin isotherm is derived using the Mie-Griineisen 
equation and the shock Hugoniot as a reference curve. An analytic solution to this 
equation is obtained in the form of power series. The arbitrary constant is expressed in 
terms of the specimen's parameters at ambient conditions. Based on this solution, the 
complete, the caloric, the thermal and the incomplete EOS are obtained. 

The proposed EOS are valid in the range of pressures and temperatures where the 
behaviour of the solid state is determined by the lattice thermal vibrations. At very 
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low and high temperatures the conduction electrons contribution must be taken into 
account. At higher temperatures the anharmonicity effects should be included as well. 

The complete EOS is derived mainly within the framework of thermodynamics 
and it cannot predict phase transitions. Should a phase transition (or melting) be 
registered experimentally, the high-pressure phase (the melt) might be described using 
the technique, given in M. 

The obtained incomplete EOS derives from a thermodynamically consistent 
complete EOS and, therefore, it might serve as a reliable physical constitutive relation, 
required to elucidate how the characteristics of the material affect the propagation of 
shock waves. This is one of the aims of the present work. 

Also, an equation for the shock temperature along the Hugoniot curve is derived. 

Without the use of the algebraic programming system REDUCE || many of these 
expressions could not have been obtained easily. 

In a forthcoming paper the proposed EOS will be applied to some solids for which 
the necessary shock-compression data are available. 
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